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Abstract 

| \ New equations are derived for meson and baryon Green's functions in a system- 

' atic method using the QCD Lagrangian and large N c approximation as a starting 

point. The equations are shown to contain both confinement dynamics and chiral 
^ , symmetry breaking. Linear Regge trajectories are demonstrated for the solutions. 

<n : 

i . 1. Introduction 

G\ ■ 
G\ ; 

The QCD dynamics of qq and 2>q systems is governed by two basic phenomena: confine- 
Q_i! ment and chiral symmetry breaking (CSB), which should be treated in a fully relativisti- 

cally covariant way. Confinement is usually introduced for static quarks via the area law 
of the Wilson loop [1] or equivalently through the field correlators in the Field Correlator 
. . i Method (FCM) [2,3]. 

.£h ! For spinless quarks, or neglecting spin-dependent mass corrections, one can envisage 

a self-consistent method which treats confinement as the area law also for light quarks in 
a relativistically covariant way. Such method was introduced originally in [4] for mesons, 
in [5] for baryons, and in [6] for heavy-light mesons, and later on in [7] an improvement 
of the method was done taking into account dynamical degrees of freedom of the QCD 
string, which naturally appears due to the area law. 

As a result Regge trajectories have been found in [7] with the correct string slope 

(27TCT)- 1 . 

Spin corrections have been considered in [8] for heavy mesons and in [6] for heavy-light 
ones, while baryon Regge trajectories have been found in [5], for review see [9]. 

In all cases the basic formalism was FCM and the Feynman-Schwinger (or world-line) 
path integral representation [3,10,11] which is well suited for relativistic quarks when spin 
is considered perturbation. 

The main difficulty which was always present in the method, was this perturbative 
treatment of spin degrees of freedom (which is incorrect, e.g., for pion) and absence of 
spontaneous CSB effects in general [12]. 

Recently a new type of formalism was suggested to treat simultaneously confinement 
and CSB and a new nonlinear equation was derived for a light quark in the field of heavy 
antiquark [13]. 
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This equation derived directly from QCD Lagrangian was found to produce linear 
confinement and CSB for the light quark [13,14] and the explicit form of the effective 
quark mass operator M(x, y) was defined obeying both these properties. 

Since the method of [13] is quite general and allows to treat also multiquark systems, 
we apply it here to the qq and 3q systems, and find dynamical equations for them, which 
contain confinement and CSB. To make our equations treatable, we systematically exploit 
the large N c limit, and mostly confine ourselves to the simplest field correlators - the so- 
called Gaussian approximation; it was shown in [13] that the sum over all correlators does 
not change the qualitative results, however the kernel of equations becomes much more 
complicated. 

The paper is organized as follows. 

In Section 2 the general effective quark Lagrangian from the standard QCD Lagrangian 
is obtained by integrating out gluonic degrees of freedom, and the nonlinear equation for 
the single quark propagator S (attached to the string in a gauge-invariant way) is derived. 

Section 3 is devoted to the qq Green's function, which can be expressed as an integral of 
product of S for quark and antiquark. Resulting equation is studied both in the differential 
form and in the integral form where again the Feynman-Schwinger representation leads 
to the effective Hamiltonian for the qq system, discussed in Appendix. 

A similar procedure is accomplished for the baryon Green's function in Section 4. 
Inclusion of chiral degrees of freedom is given in Section 5, where an effective chiral 
Lagrangian is derived for bosons. Discussion of results and comparison to other methods 
is done in Conclusion. 

2. Effective quark Lagrangian 

As was discussed in the previous section, one can obtain effective quark Lagrangian by 
averaging over background gluonic fields. We shall repeat this procedure following [13] 
now paying special attention to the dependence on the contour in the definition of contour 
gauge, and introducing the operation of averaging over contour manifold. The QCD 
partition function for quarks and gluons can be written as 

Z = J DAD^D^ + exp[L + L x + L int ] (1) 

where we are using Euclidean metric and define 

Lo = -\j d'x(F^) 2 , (2) 

L 1 = -iJ f ip + (x)0 + m f ) f ij;(x)d 4 x, (3) 

= J f ^ + (x)gA(x) ^{x)d A x. (4) 

Here and in what follows ^ip aa denotes quark operator with flavour /, color a and 
bispinor index a. 

To express A^x) through F^ u one can use the generalized Fock-Schwinger gauge [15] 
with the contour C(x) from the point x to Xo, which can also lie at infinity, 
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Now one can integrate out gluonic field A^x), and introduce an arbitrary integration 
over the set of contours C(x) with the weight D K (C), since Z is gauge invariant and does 
not depend on contours C(x). One obtains 

Z = J £>K(C)£ty£ty+exp{Li + L cS } (6) 

where the effective quark Lagrangian L c s is defined as 

expLcfr = (exp J f i/; + A f ^d 4 x) A - (7) 

Using cluster expansion L e g can be written as an infinite sum containing averages ((A) k )A- 
At this point one can exploit the Gaussian approximation, neglecting all correlators {(A) k ) 
for degrees higher than k = 2. Numerical accuracy of this approximation was discussed 
in [9,16] and tested in [17]. One expects that for static quarks corrections to Gaussian 
approximation amount to less than 10%. 

The resulting effective Lagrangian is quart ic in ip, 

= ^ / d A xd A y^t a {x) f ^(x) ^ft(v) ^ aS (y)J a ^s(x, y) + 0(^ 6 ), (8) 

y) = {i^) a p{iv) 1 5J^{x,y) (9) 

and J^y is expressed as 

J^v{x,y) = 9 2 J c ]^T du e J c -^ du e'^( F eM F e'u>{v)) . (10) 

L cS (8) is written in the contour gauge [15]. 

It can be identically rewritten in the gauge-invariant form if one substitutes par- 
allel transporters xq), xq) (identically equal to unity in this gauge) into (8) 
and (10), multiplying each ip(x) and ip(y) respectively and in (10) replacing F(u) by 
u)F(u)$(u, Xo) and similarly for F(v). 

After that L c s becomes gauge-invariant, but in general contour-dependent, if one 
keeps only the quartic term (8), and neglects all higher terms. A similar problem occurs 
in the cluster expansion of Wilson loop, when one keeps only lowest correlators, leading 
to the (erroneous) surface dependence of the result. 

Situation here is the same as with a sum of QCD perturbation series, which depends 
on the normalization mass \x for any finite number of terms in the series. This unphysical 
dependence is usually treated by fixing /i at some physically reasonable value /x (°f the 
order of the inverse size of the system. 

The integration over contours Dk(C) in (6) resolves this difficulty in a similar way. 
Namely, the partition function Z formally does not depend on contours (since it is inte- 
grated over a set of contours) but depends on the weight Dk(C) and we choose this weight 
in such a way, that the contours would generate the string of minimal length between q 
and q. Thus the physical choice of the contour corresponds to the minimization of the 
meson (baryon) mass over the class of strings, in the same way as the choice of \x = fiQ 
corresponds to the minimization of the dropped higher perturbative terms. 
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As a practical outcome, we shall keep the integral Dk{C) till the end and finally use 
it to minimize the string between the quarks. 

Till this point we have made only one approximation -neglected all field correlators 
except the Gaussian one. Now one must do another approximation -assume large N c 
expansion and keep the lowest term. As was shown in [13] this enables one to replace 
in (8) the colorless product ^ipb( x ) 9 ipb(y) = ty( x )®( x i x o)®( x Oiy) 9 4 )+ (.y)) by the 
quark Green's function 

'iM*) 9 ^(y) - 5 fg N c S^(x,y) (11) 

and assumes the form 

= -ij d'xd'y ty+ (x) ^M aS (x, y) ^ aS (y) (12) 
where the quark mass operator is 

f M aS (x,y) = -J iUI/ (x,y)(7 jU f S(x, y)^ v ) aS . (13) 
From (6) it is evident that *S satisfies to equation 

(—id x - irrif) f S(x, y)-ij f M{x, z)d 4 z f S(z, y) = 5 (4 \x - y). (14) 

Equations (13), (14) have been first derived in [13]. From (6) and (12) one can realize 
that at large N c the qq and 3q dynamics is expressed through the quark mass operator 
(13), which should contain both confinement and CSB. 

Indeed, analysis done in [13,14] reveals that confinement is present in the long-distance 
form of M(x,y), when both distances |x|, |y| of light quark from heavy antiquark (placed 
at x = 0) are large. 

We shall do now several simplifying assumptions, to clarify the structure of M(x,y). 
First of all we take the class of contours C going from any point x = (x4,x) to the 
point (x4, 0) and then to (— oo) along the For this class the corresponding gauge 

was studied in [18]. Secondly, we take the dominant part of in (13), namely J44, 
which is proportional to the correlator of color-electric fields, yielding linear confining 
interaction, and neglect other components .1^ J i4 , Ja,i = 1,2,3, containing magnetic 
fields and yielding momentum dependent corrections. (It is easy to take into account 
these contributions in a more detailed analysis). 

The correlator (FF) in (10) can be expressed through the scalar correlator D(x), 
defined as [2], 

-^{F af3 {u)${u,v)F lS (v)$(v,u)) = D(u - v)(5 ai 5p S - <WW> + O(A) (15) 

where the correlator D\ not contributing to confinement is neglected. As a result one has 
for M [19] 

f M Cxi (x,y)= f M^I+ f M^a t + / #' 74 + f M^^. (16) 
Here we have defined 

»< = (? '1 < 17 > 
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The dominant part of M, f is linearly growing at large |x|, |y| and in the most simple 
case of Gaussian form of D(x), can be written as 



fM<®(x,y) = — -=e -^- a \^\~5^(x-y) (18) 

where T g is the gluon corretion length, and 5 is a smeared 5-function, which can be 
represented as [19] 

5( 3 )(x-y)^e^(-^^)(^=) 3 , b~2T g . (19) 

Here T g is the gluon correlation length, which enters D(u) as D(u) = D(0)exp(— ^-). 

We are now in the position to derive qq,3q Green's function, which will be done in the 
next section. 

3. Equations for qq Green's function 

We start with gauge-invariant definitions of initial (^m) and final (^&n) qq states, 

V- m (x,x) = H + (x)r'f ih <Z>(x,x) H(x); 

*&n(y,y) = f ^ + (y)rf; f -^(y,y) Hiv)- (20) 

The qq Green's function is expressed as 

G(xx\yy) = (V- m (x, x)^ n (y, y)) = 

= 1 J DK(C)DijDij + exp{L + L eff }* in *+. (21) 

One can do integrals over D^Dijj + and assuming as before the large N c limit, one can 
neglect the determinant term, yielding 

G(xx\yy) = 1 1 Dk(C)MS§(x, y)r££s£(y, x)r% fi }- 

-tr(s£(x,x)r™)tv(si(y,y)r + ^)} (22) 
where we have omitted for simplicity parallel transporters and used notation 

S f c (x, y) = (-id - im f - iMf(x', y'))^. (23) 

At this point we assume that minimization over contours Dk(C) in (22) selects one specific 
contour C{u) for each u (and we find this C(u) in what follows), and choose for simplicity 
nonsinglet flavour channel. Then the second term in the curly brackets in (22) disappears 
and we have equation for G 

(-id x - imi - iM Cl (x,y))(-idx - im 2 - iM C2 (x,y))G(y,y\z, z) = 

x - z)5^\x - z) (24) 
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where we have omitted integrals over d 4 yd 4 y and unit operators 1 = S^ 4 \x — y),l = 
S^(x — y) as factors of (— id x — imi) and (— id^ — im 2 ), respectively. Next steps in 
treating (24) are standard for the Bethe-Salpeter formalism [20], the only difference is 
the contour dependence of mass operators Mq, which we specify as follows. One can 
choose the class of contours which pass from a point u = (u^, u) to (u^, x ) and then to 
(— oo,x ), i.e. parallel to the cm. motion of the qq system. Exact value of x will be 
found by minimization of the interaction kernel. 

Defining as in [20] the cm. coordinate X M and relative coordinate r M . 

X M = ax^ + (1 - a)xn, r lk = x ik - (25) 

where a is an arbitrary parameter and corresponding momenta 

* - \h p » - \k (26) 

one can also fix a to have correspondence with nonrelativistic limit to be a = 
mi (mi + m 2 ) _1 . 

Since Mq is invariant with respect to time shifts, one can introduce the total mass of 
the qq system E and relative energy e and write in the cm. system an equation for the 
qq wave function 

1 ; " : £-(mi/3i + p-ai)+e-/3itf(x-Xo)] 



mi + m 2 



x [—^ — E - {m 2 {3 2 - p • a 2 ) - e - &l/(x - x )]^ = (27) 
mi + m 2 

where C/(x — x ) is the local limit (T g — > 0) of the mass operator M c . For vanishing 
angular momentum one has 

C/(x-x ) =(t|x-Xo|. (28) 

The further analysis and solution of (27) involves as in [21] introduction of positive 
(+) and negative (-) energy projection operators and corresponding 4 wave functions 
\I/aa'> \ A' = ±, for which a system of equations is written similarly to [21]. 

We shall not follow this route here, however, leaving it for subsequent paper and 
instead write in the Appendix the Hamiltonian approach to the same problem, which will 
allow us to calculate spectrum in an approximate way, and discuss it in the Conclusions. 

There is still another the 3d form of (27) which obtains in the same way as in [21]. 
Indeed, writing the time Fourier trasform of the Green's function in the form (the sign 
()c denotes the integral Dk{C)) 

^ (E — E 2 — H\)(E 2 — H 2 )^ C ^ ^ 

and integrating over dE 2 one arrives as in [21] to the 3d Green's function 

G(S,r,rO = ( 1 ) c (30) 
hi — rl\ — n 2 

where 

H 1 = rriiPi + p • cki + /3iC/(x - x ), 
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H 2 = m 2 /3 2 -p-a 2 + /3 2 [/(x-Xo). (31) 

Minimizing interaction over x , one arrives at x = ^rp^, and at the equation for the 
wave function. 

{mi/3i + m 2 (3 2 + p(ai - a 2 ) + + /W-^)}*(r) = £tt(r). (32) 

In the nonrelativistic limit this reduces to the usual nonrelativistic quark model, 

^P 2 , ^ t? m x m 2 

(- h ar)ip n = e n tp n , e n = E n - m 1 - m 2 , m = ■ (33) 

2m mi + m 2 



4. Equations for the baryon Green's function 

Equations for the 3q system can be written in the same way as for the qq system. We 
again shall assume large N c limit in the sense, that 1/N C corrections from qq pairs to the 
quark Green's function and the effective mass can be neglected, but shall write explicit 
expressions for N c = 3. 

The initial and final field operators are 

* in (x,y,z) = e abc T a ^^ aa (x,C(x))^(y,C(y))^ n (z,C(z)) (34) 

with the notations: a, b, c, are color indices, a, (3, 7 are Lorentz bispinor indices and trans- 
ported quark operators are 

Vw(z, C{x)) = ($ c (x, x)ij) a (x)) a (35) 

and the contour C (x) in $<7 can be arbitrary , but it is convenient to choose it in the same 
class of contours that is used in Dk(C) and in the generalized Fock-Schwinger gauge [15]. 
pa/37 i s the Lorentz spinor tensor securing proper baryon quantum numbers. One can also 
choose other operators, but it does not influence the resulting equations. In (34) we have 
omitted flavour indices in T and if)(x, C), to be easily restored in final expressions. 

Using now the effective Lagrangian (12) valid at large iV c , we obtain for the 3q Green's 
function. 

G^\x,y,z\x',y',z') = 
^ J Dk(C)D^D^ &q (x\ y', zf)*£(x, y, z)exp(L l + L cS ). (36) 
Integrating out quark degrees of freedom and neglecting determinant at large N c one has 

G {3q) = J DK(C)(er)(eT'){S(x, x')S(y, y')S(z, z') + perm} (37) 

where color and bispinor indices are suppressed for simplicity together with parallel trans- 
porters in initial and final states. 

One can also define unprojected (without T, V) 3q Green's function Gffl with 3 initial 
and 3 final bispinor indices instead of projected by T, V quantum numbers of baryon. 
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Assuming that minimization over contours Dk(C) reduces to the single choice of the 
contours (the single string junction trajectory minimizing the mass of baryon), one can 
write equation for Gffl: 

(-id x - im x - iM x ){-idy - im 2 - iM 2 )(-id z - im 3 - iM 3 )G { ^ ] = 

x — x')5^\y — y')8^\z — z') (38) 

and e.g. M ± G = J M(x,u)G(u, x')d 4 u. One can simplify the form (37) for G {3q) taking 
into account that M(x,x') actually does not depend on ^£±5^ an d hence the interaction 
kernel of G^ does not depend on relative energies, as in [21]. Similarly to [20,21] one 
can introduce Fourier tranform of G^ in time components and take into account energy 
conservation E = E\ + E 2 + £"3. One obtains 



G^(E,E 2 ,E 3 )~ f DK{C){eY){eV) 



1 (39) 



(E — E 2 — E 3 — H 1 )(E 2 — H 2 ){E 3 — H 3 ) 
where notations are used 

H t = mS {i) + p (i) a« + /3«M(r« - r^) (40) 

and we have taken in M(x, x') the limit of small T g and set of contours in Dk(C) passing 
from the point to some (arbitrary) point r( ). 

As in [21] one can now integrate over E 2 , E 3 to obtain finally 

r^) ~ / D < C )W^') {E _ Hi l _ H2 _ Hz y ( 41 ) 
From (41) one obtains equation for the 3q wave function similar to that of qq system, 

(H, + H 2 + H 3 - E)ii(r u r 2 , r 3 ) = (42) 

where r*- -* is to be taken at the Torricelli point. 

In the nonrelativistic approximation rrii ^> yfa one has 



E 



IP I 1 a \ r (i) _ r (0)| 

2m. 



# e = E-J2rm (43) 



5. Inclusion of chiral degrees of freedom 

Analysis of the qq equations, e.g. (24) and (32), reveals that solutions do not contain 
Goldstone modes in the chiral limit (mi,m 2 — > 0). The reason is that the equations have 
been obtained in the limit of large N c , while chiral corrections appear in the subleading 
order, and the coupling constant of pion to the quark is 0(N~ 1 / 2 ). Therefore chiral 
correction 0(1/N C ) to the second term on the r.h.s. of (22) yields contribution of the 
same order 0(N C ) as the first term, since number of color traces in these terms are 
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different. A similar situation occurs in the vacuum made of instantons surrounded by 
confining background considered in [22]. There the massive p-type pole in the pionic 
channel produced by the first term in (22) is exactly cancelled by the second term, while 
the Goldstone pole appears in the chiral correction. 

To proceed one can bosonize the quartic Lagrangian (8) in the standard (but nonlocal) 
way [23] 

Z = J DK^D^+Dioexpi^ + L u + L^} (44) 
where bosonic colorless field uo = ur^(x, y) enters 

in L^, as 

1 r 

L u = -— J d 4 xd 4 yu f a 9 3 ( y x,y)J a(3 .^( y x,y)u f 1 9 5 ( y x,y) (45) 

= J d A xd 4 y[uj f a 9 p {x,y)J a p. lS {x,y) f rf(x) 9 if> s (y) + h.c.]. (46) 

Integrating out quark degrees of freedom one obtains the effective chiral Lagrangian, 
Z ~ Je L ^DLuDK(C) 

L c h = -L u + Trln(i<9 + im + %M + A) (47) 

where M is the leading 0(N®) term (13), while A is the chiral correction. 

Note that in the limit T g — ► 0, the quark mass operator M becomes local (~ S^(x — 
y)) and therefore one can consider in A and the usual limit of local chiral field. 
Parametrizing it as in [24], one has 

A u = ?(e i7rT75 - 1)M. (48) 

Proceeding as in [22] one can expand L ch (7r)in 7Tj and obtain 

L ch (n) = J n a {k)N{k,k')Ti a {k')dkdk' (49) 

where N(0, 0) = yielding massless Nambu-Goldstone pions. Similarly to [22] also 
here the massive pole cancels while the Nambu-Goldstone pole survives. For a detailed 
discussion the reader is referred to [22] and a subsequent publication. 

6. Conclusions 

We have obtained above equations for mesons and baryons (27), (32) and (38), (42) 
respectively which contain all dynamics in the form of nonlocal mass Kernels M(x,y). 
The latter are solutions of equations (13), (14) obtained earlier in [13] and describing 
motion of one quark in the field of heavy antiquark. This reduction of qq and 3d problems 
to a simpler problem of one quark in the field of static source is possible in the lowest 
order of 1/N C expansion. 

The form of solution for M(x,y) was found quasiclassically in [13,14] and also numer- 
ically in [19] in the thin string limit T g — > 0) and for nonrotating string (L = 0) reduces 
to a simple linear potential at large distances. 
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In the same limit the obtained equations produce reasonable spectrum, see, e.g., equa- 
tion (A. 12) for the qq system, which yields results roughly coinciding with the WKB 
spectrum for spinless Salpeter equation in a case of zero current quark mass [25]. 

Thus confining dynamics is correctly reproduced by the cited equations. At the same 
time the exploited form of M(x, y) is applicable at zero angular momentum and does not 
describe the rotating string, which can be seen in the wrong Regge slope in (A. 12). Chiral 
degrees of freedom which can be taken into account in the same formalism need additional 
step, which was shortly discussed in Section 5. Here one should take into account the 4q 
form of the effective Lagrangian (8) and do a general (nonlocal) bosonisation procedure, 
which generates a new operator M(x,y), containing bosonic fields, in addition to the 
scalar-isoscalar component of M, dominant at large N c , which represents the scalar string 
(and hence chiral symmetry breaking). 

The paper has left aside many interesting questions, in particular the exact form of 
equations for pion Green's function, the coexistence of confining and chiral effects in the 
dynamics of mesons and baryons, the definition of the constituent quark mass and hadron 
magnetic moments etc. 

All these points are planned for future publications. 

The author is grateful to Yu.S. Kalashnikova, A.V.Nefediev and V.I.Shevchenko for 
useful discussions. 

The financial support of RFFI through the grants 97-02-16404, 97-02-17491 and 96- 
15-96740 is gratefully acknowledged. 
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Appendix 

Effective Hamiltonian for the qq system 

Our purpose here is to use the Feynman-Schwinger representation (FSR) [10,11] to express 
the inverse operator in (23). To this end we represent Sq as 

S f c (x,y) =i(d + m + M)- 1 = i(-d + m + M)(-d 2 + 5 + M 2 y l (A.l) 

where 5 = dM — Md, and with the notation M 2 = M 2 + 5 one has using FSR 

poo 

S f c (x, y) = i(-B + m + M) / dse~ s{ - d2+M2) = 

Jo 

i(—d + m + M) ds(Dz) xy e- K -^ dT ^ {zr) . (A.2) 
Proceeding now in the same way, as in [7] one obtains 

G(x,x\y,y) = ±j Dk{C)Y + {-8 + m h + M fl )F(-d + m fl + M fl ) 

x J DRDrDfi 1 Dfi 2 e~ A (A.3) 

where we have used notations 

R a = !^L±!^L i r a = z a (t)-z a (t),t = z A , (A.4) 

dz 4 dz4 
2/ii = — ,2// 2 = 

dr dr 

The Euclidean action A can be written as 

A= dtl + 2 — + ^-^R 2 + 

Jo \2/ii(t) 2fi 2 (t) 2 2 

M?(z(t)-x ) M s (z(t)-x ) \ 
+ 2^ + 2^ }■ (A ' 5) 

with the notations \x = , /^+ = A*i + /i 2 - Integration over Dk(C) in (A.3) can be 
reduced to the minimization with respect to the parameter (3, defining x (/3) = z/3 + 
z(l -P),0 < (3<l. 

At large r one can use the asymptotic form of M, given in (18), then neglecting for 
the moment the term S which is constant at large r one has 

tt2 M 2 (z-x°) M 2 (z-x°) 22/ (l-/9) 2 (3 2 , IK x 

2//i 2/i 2 v 2/ix 2^ 2 ; v ; 

Finding the stationary point of U 2 in (3 at /3q = — , one finally obtains the minimized 
U 2 , denoted as < U 2 > 



a 2 r 2 



< U 2 >= (A.7) 

2/iH 
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Integrating over DR^ with boundary conditions R±(T) = T and R±(0) = 0, one has for A 
The corresponding Hamiltonian is 

2 2 2 2 2 

mi mi a + crr z p z . . _ 

#= — + — + — + + — . (A.9) 

2/i! 2/i 2 2 2/1+ 2^ v ; 

One can vary (A.9) over //i,//2 with fixed /t to obtain for equal masses m± — m 2 = m 
and finally varying over \x one has 



H = + m 2 + a 2 r 2 (A.ll) 

which yields the mass of the qq system 

M 2 = 8a(2n r + L + |) + 4m 2 . (A. 12) 

It is remakable that the spectrum (A. 12) is the same (modulo factor of 4) as the WKB 
spectrum of Dirac equation with the linear confining potential [13]. 

However the Regge slope of the spectrum (A. 12) is (8o") _1 , which signifies that the 
rotation of the string connecting q and q is not properly taken into account. Here and 
above we confine ourselves to the zero angular momentum, leaving the problem of rotating 
string to subsequent publication. 
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